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Abstract. If Kermit's What-Happens-Next-Machine had functioned, you would not 
have seen much, because it would have gone too quickly. In this article it is shown 
that putting up and solving the equations of motion of a seemingly simple mechanical 
apparatus presents a challenging problem. The simulation can, however, be quite 
instructive and also entertaining. 



1. Introduction 




Fig.l Scheme of the What-Happens-Next- Fig. 2 Geometry of the seesaw and the box's lid. 
Maschine. Inset: Kermit the Frog in the scene. Definition of positions, distances and angles. 

In the children's television series Sesame Street Kermit the Frog, a muppet created 
by Jim Henson, presented his What-Happens-Next-Machine. Although he proclaimed 
having constructed it to avoid walking from here to there to turn on his radio, the 
purpose of the scene, of course, lies in conveying to the spectator the concept of a 
chain reaction. More spectacular chain reactions are now often demonstrated in science 
museums and TV shows |2]. The What-Happens-Next-Machine, however, is a - if not 
the - historical, famous or notorious example. The setup is meant to function as follows: 
A string being cut, a sandbag hung from a pulley drops onto a seesaw. To the other end 
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of the seesaw a rope is tied which hfts the hd of a box containing a balloon, preferably 
filled with helium. From the balloon a string had already been trailing out of the box 
and is attached to the switch of a radio. So as the lid opens and the balloon rises, at 
the point the strings is pulled straight up the radio is turned on. But, fitting the typical 
irony of Sesame Street, even though Kermit tries to save whatever possible, nearly every 
step - except for the balloon floating up (with some delay) - is malfunctioning pj. Thus, 
albeit for the good entertainment the flawed performance provides, unfortunately the 
spectator does not get to see the before described and intended chain reaction. Without 
malfunctions, how fast would it go? 




Fig. 3 What-Happens-Next-Machine built by 
students in Konstanz 2009. Slightly difFent from 
the original and theoretically discussed version, 
the rope is attached at the middle of the box's 
lid and the seesaw's little edge will prevent the 
sandbag from sliding. The inset shows a rocker 
switch a straightening ballon string can actually 
operate.) 




Fig. 4 Here the sandbag has fallen, turned over 
the seesaw and opened the lid. In this foto 
taken during setup there is no balloon yet. The 
seesaw need not pull the lid up to the vertical. 
The balloon can be smaller than the box and 
come out sideways. Anyway, if the lid is pulled 
strongly enough at the beginning, it has enough 
momentum to move on up to vertical. 



While it can also be a nice task for students to build and run the machine, this 
article focusses on calculating the motion. Referring to every-day life phenomena or 
using toys |1] can make physics courses and exercises more appealing. The difficulty 
often is that one physical concept in its pure form is insufficient to explain the behaviour. 
A calculation of the What-Happens-Next-Machine as well has to begin by making 
appropriate idealizations. It can, however, provide a challenging exercise for second- 
year students of physics. In analytical mechanics textbook examples may fall short of 
actually conveying the virtue of the Langrange formalism, because the posed problems 
of calculating something's motion can just as well be solved by Newtonian concepts, 
evaluating forces and using energy and momentum conservation. In the setup considered 
here in the connection of the seesaw to the also rotating lid by the rope an unexpected 
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complexity is concealed. Finally, numerics will be used making the theoretical problem 
also a programming exercise. 

2. Lagrangian 

The setup is drawn in Fig.l and a close-up with labels is given in Fig. 2. We need: 
fixed quantities: 

g : constant of gravitation, 9.81 m s"^ 

ms : mass of the sandbag (considered a point mass) 

rn : length of the right half of the seesaw 

: length of the rope (straightened) 
ri : half edge of the box, and thus also of the lid (assumed square) 
mi : mass of the box's lid 

II : moment of inertia of the box's lid with respect to its hinge; II = w-l ■ (2rL)^/3 

tb : radius of the balloon (assumed a sphere); we set rs = 

PAir '■ density of air, 1.29 kg m^'^ pHe- density of helium, 0.18 kg m~^ 

variables: 

xs : horizontal position of the sandbag counted from the middle of the seesaw 

hs : vertical position (height) of the sandbag also counted from the middle of the seesaw 

rs : distance of the sandbag from the middle of the seesaw when the two are in contact 

ifs : tilt angle of the seesaw, zero=horizontal, pos.=left side is down, neg.=left side is up 

cpL : opening angle of the lid, counted positive when opened from horizontal 

h B : height of the balloon (center) 

The following energies are involved: 

component potential kinetic 

sandbag ms g hs ^ (i;| + /i|) 

lid rriL g sin Lp^ ^ <p1 

balloon ^ir r| (pHe - PAir) 9 hB 5 |vr r| pHe h\ 

hs=—rssm(fs if the sandbag is on the seesaw. The kinetic energy of the sandbag 
can be rewritten as ^(r| + rgip'g) if the sandbag and the seesaw are moving together. 
The potential energy of the balloon accounts for buoyancy. The seesaw is regarded as 
massless (that is negligible mass as compared to the sandbag). The box is fixed to the 
ground, only the lid moves. The rope is assumed massless and the mass of the balloon's 
rubber skin is neglected, too [5]. There is no friction in the hinges of the seesaw and 
the lid and as well no sliding friction is considered should the sandbag slide down the 
seesaw. The Lagrangian [6] (for the sandbag in contact with the seesaw) is therefore 
put up as: 

>C = — (r5 + rg^s) + "^sgrs sin Lps + y 0l - "^Lfl"^L sin v^l (1) 

14 -4 

+ ^^T^r^pHehl - -vrr|(pHe - pAir)ghB 

In the following we shall, however, consider C without the balloon's contributions. To 
neglect the pressure the balloon due to buoyancy exerts from below onto the lid, seems a 
realistic assumption. The balloon will therefore only follow the motion of the lid if it can 
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rise fast enough. We further assume that the balloon does not bounce off downwards 
when it hits the lid, but that its motion is completely damped out in this case. Letting 
aside the balloon, following the Lagrange formalism we deduce: 

.dC d dC . dC d dC . dC d dC 

ors dtors ocps dtocps uLpi dtOLpi 

{msrsVs + '^sQ sin v^s - msfs)5rs + {msgrs cos ips - '^msrsrs^s - 
msrl(ps)5^s " {migrL cos ipi + h(pL)5^L = (2) 
t is time, rs is an independent variable, whence we get 

f 5 = g sin ips + rs 0l (3) 

as a first differential equation or equation of motion. With the rope straight, (ps and (pi 
are not independent. We choose to take (p^ as a function of (ps- The geometric relation 



is given in an attachment. Then S(pL = :^^Vs, <Pl = ^^05 and (pi = -^^s + ^Vs- 
The derivatives ^ and '-^r^ are evaluated numerically in an actual calculation. Thus 

dtps dipj, J 

a second equation of motion follows from (2): 

msgrs cos ips - 2msrsfs^s - (rriLgrL cos (pi + h ^ 0|) ^ 

msri + lLi^y 

The functionality ipiiips) is a constraint in the sense of the Lagrange formalism. It 
takes the form of an equality if the rope is straight. A non-holonomous constraint, 
nevertheless, is also always present if the rope is loose, namely that the right endpoint 
of the seesaw and the right edge of the lid cannot be more than a distance of Ir 
apart. Even if here we do not handle the balloon as part of one complete Lagrangian, 
there is an also non-holonomous constraint on its motion. It cannot exceed a certain 
height determined by the lid, whereas it may well stay below that limit. Geometric 
considerations concerning this point are given in an attachment. Furthermore, not 
including the balloon pressing against the lid from underneath into the equations of 
motion, care must be taken not to let the lid drop below its initial position on top 
of the balloon. The integration of (3) and (4) is performed numerically, in the most 
simple manner in the attached program. More advanced methods like, for example, the 
Runge-Kutta scheme, could of course be employed. 



3. Impact 

Calculating the free fall of the sandbag as well as the time and position of its hitting the 
seesaw is trivial. At the beginning the seesaw can be at rest at any angle the length of 
the rope allows. As long as or if the rope is not straight our massless seesaw will just give 
way to the free falling sandbag. However, the instant the rope straightens or is already 
straight when the sandbag hits the seesaw, because of suddenly the weight of the lid 
pulling there will be a discontinuity in the velocities, which must be determined apart 
from equations of motion like (3) and (4). This turned out to be the most difficult, and 
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to my opinion still somewhat ambiguous aspect of the problem. For reasons explained 
later the more general case is investigated that the lid can already be moving when the 
rope straightens or the sandbag hits the seesaw. Furthermore can the sandbag also have 
a horizontal velocity in its free fall. As there is nothing to stop it, the velocity of the 
sandbag along the seesaw v\\ is conserved and assigned to rs at impact. 

rs = -xscosips - hs simps (5) 

But what about the initial angular velocity of the seesaw immediately after impact? We 
shall assume an inelastic collision insofar as the sandbag never bounces off upwards from 
the seesaw. Immediately after impact the seesaw and the sandbag move together (later 
it is checked whether the vertical velocity from the turning and centrifugal motion 
would exceed free fall and eventually the sandbag is allowed to loose touch with the 
seesaw). Both kinetic energy and angular momentum cannot be conserved in the 
collision. Drawing the analogy to the one-dimensional inelastic collision where the two 
bodies stick together afterwards, momentum is conserved, but part of the kinetic energy 
goes into deformation and heat, the logical decision in our case is to require conservation 
of angular momentum. However, as the sandbag on the seesaw and the lid do not turn 
around the same center, the definition of angular momentum is not necessarily clear here. 
It further only imports at the singular moment of impact. (The motion afterwards of 
all parts together as calculated from the Lagrangian obeys energy conservation.) 

Our treatment of the impact is therefore based on the following very elementary 
considerations [U The instant the sandbag hits and/or the rope straightens, the pulling 
lid acts like a hypothetical mass on the seesaw at distance rjj from the center. The 
angular momentum the lid has with respect to its hinge can instead be taken as that of 
one third of the mass sitting at its end. Therefore decompose Ii'^LMfore into a lever arm 
2rL and momentum perpendicular to it p^^i = J^^^^^^^I^_ However, now in order to find 
the part having an effect at the other end of the rope tied to the seesaw, the momentum 
has to viewed as made up of two in general non-perpendicular components along the lid 
and along the rope, the latter being pr^l = -^^^f^- How to evaluate the angle a is also 
given in an attachment. At the seesaw this momentum translated along the rope with 
lever arm th makes up an angular momentum ^'"'^^'''^^"'''^'^"^^^^'^'^^ '^^^ (Fig. 5). After 

I Nevertheless, conserving the velocity of the sandbag along the seesaw and for the rotations angular 
momentum in the form msr^ips+lL^L does not produce spurious results (rs and ips before the collision 
just describe the sandbag at the end of its free fall in polar coordinates, (pL before is known, afterwards 
'PL = Vl{'^s,Vs), thus only tpg needs to be determined). Neither does conserving kinetic energy of 
rotational motion ^^r'gip^g + produce any qualitatively or obviously impossible behaviour (with 

in general the rope at oblique angles to the seesaw and the lid, forces perpendicular to one can well 
be along the other, and thus conservation of the sum of angular momenta of each around their hinge 
need not necessarily be expected, which would leave energy as a quantity eventually to be conserved). 
Attempting to evaluate velocities after the collision by normal momentum conservation together with 
the constraint that afterwards the sandbag on the seesaw and the lid at least momentarily move together 
does lead to contradictions, even if the angles '^s{t) and V'i(i) do not show them. There is more kinetic 
energy after the collision than before. And this does produce cases where the sandbag bounces off the 
seesaw with an upward velocity component. 
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the collision the change of (fi is determined by that oi (ps, nevertheless, the "transfer 
coefficient" for angular momentum through the rope remains the factor sm(a+yg-yL) ^ 
With v± = xs sin ips — hs cos ips the velocity of the sandbag perpendicular to the seesaw- 
before, to solve for ips immediately after impact we have the following equation: 

h<fL,beforerH sm{a + p>s - V^l) _ 



2ri sin a 



{msrl 



IlTh sin(a + (ps - 'Pl^ — 
2rL sin a 



dips 



(6) 



The first time the sandbag hits when it comes from vertical free fall and the lid is still 
at rest, we have xs = and ^pLfiefore = at impact. Of course, the argument could 
have been led the other way round, first replacing nis at by some mass turning at 
vh on the seesaw, translating a momentum component down the rope to the lid and 
establishing conservation of angular momentum around the hinge of the lid. The results 
for (fs and (fL after the collision would have been the same. 
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Fig. 5 Determining angles and evaluating the Fig. 6 Gravitational and acceleration forces and 
momentum component of the lid along the rope. tearing forces on the rope. 



4. Torques, forces and energy 

In this section it shall be sought whether the equation of motion could have been deduced 
by more elementary means than the Lagrange formalism. With the sandbag on the 
seesaw the sum of kinetic and potential energies is 

/ 



E 



— (r + rgLps) - msgr-s smLps + y^l + mLgri sm ^pl 



(7) 



Requiring energy conservation leads to: 
dE 







dt 



msTsTs - msTsrs'Ps - i^sg^s smp>s + '^rnsrsVs'Ps 



(8) 



+ msrlips^s - msgrsips cos (ps + Il^l^l + migri^L cos ipi 

where the term msrs'f's'^'s been split into a negative and twice the positive 
contribution. Now using the knowledge that the radial motion of the sandbag will 
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follow the respective component of the gravitational force and the centrifugal force, 
futher demand that the first three terms in (8) and the rest be zero each on its own. 
Dividing by msrs or ^5, respectively, and expressing ipi and by derivatives and ips 
and (fs as before for the motion with the rope straight, this will reproduce (3) and (4). 

However, as likely if not more likely than by putting up energy conservation, one 
would approach a mechanics problem with levers and ropes by looking at forces and 
torques. The following considerations are still restricted to the motion when the sandbag 
is in touch with the seesaw and the rope is straight. As explained earlier, the first 
equation of motion (3) for is rather obviously written down. The torques determining 
the angular accelerations (fs and (pi are more difficult to find. The time derivative of the 
sandbag's angular momentum msr'gips is msr'gip s + '^^sf s'f' s'P s ■ Like in the last section 
the change in angular momentum of the lid is taken into account as seen at the other end 
of the rope. Therefore with respect to the seesaw's center this is Il'^l '^^[■""sina"'^^'' • "^^^ 
component of the gravitational force on the sandbag which leads to angular acceleration 
is msgrs cos(ps- For the torque m^'s gravitation exerts on the lid replace the the lid 
by half the mass ^ at the end of an otherwise massless lid. As force component 
perpendicular to the lid there you have coslpl, multiply by 2rL to get the torque. 
This torque translates into an effective one cos ipL ■ ■ ^" '^^[."ti'na"'^^^ other 
end of the rope up at the seesaw, which is, of course, opposed to the sandbag's. The 
seesaw and the lid are literally tied together and therefore the available torques together 
have to provide the changes in angular momentum taken together: 

th sin(a + <fs - Vl) 



2rL sin a 



mi „ th sin(a + v^s - V^l) 

msgrs cos ^ps - —g cos v?l ■ 2rL -. (9) 

2 2rL sm a 

We checked numerically that ^ = "''''2r"tin«"'^^^ • when again = ^<fs + 

^j^^s is implied, (9) becomes identical to (4). Setting up the equation of motion for 
ips from forces and torques one is likely to forget the term 2msrsTs'pSi which looks 
like a Coriolis force [7J, because (3) already seems to cover effects associated to the 
radial motion. Instead of projecting the lid's change in angular momentum into an 
acceleration of rotational motion of the seesaw, one could decompose forces that lead 
to the respective torques at the right end of the seesaw and at the end of the lid into 
parts along the seesaw or lid and along the rope and say that the tearing force on the 
rope is equal in both directions (Fig.6). The acceleration forces are against the sense of 
the angles, but an equivalent force on the right side of the seesaw to the weight of the 
sandbag on the left has to point upwards. Further Ts/th is for changing the lever arm 
of the sandbag from rs to th and sin(a + ~ ^l) and sin a are the already known 
geometrical factors from the parallelograms of forces (or momenta) to get parts along 
the rope. 

1 frriL , Il(Pl. msgrs cos ^ps - msrlips - 2msrsrs^s 

gcos(pL + ^ — ) = — , ^ (10) 



sin a 2 2rL rn sm(a + (ps — '^l 

(10) is the same as (9). 
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read parameters and initialize values 

especially set (sandbag on seesaw) to FALSE, and (rope straight) to FALSE. 



begin time loop 



check if Cpg or (p^ has reached 71/2 
i NO 



YES 




check (sandbag on seesaw) 
I FALSE ~ 



TRUE 



sandbag falls freely (keep iiorizontai velocity), lid falls freely or rests on balloon 

check if rope straightens (distance end of seesaw- end otUd) 

-* eventually set (rope straight) to TRUE and seesaw follows lid 



check if sandbag hits seesaw, if yes: set (sandbag on seesaw) to TRUE 



NO 



I YES 



if (rope straight): reinitialize velocities due to impact 



if (rope straight): 
sandbag on seesaw 
and lid move together 



if not (rope straight): 

sandbag pushes seesaw out of its way 

lid falls freely or rests on balloon 



check if sandbag's downwards velocity exceeds free fall 
if yes: set (sandbag on seesaw) and (rope straight) to FALSE 
and do not count this time step 



^ NO 



YES 



if (rope straight): check force 
direction (connection-rod check) 

-» eventually set (rope straight) to FALSE 



if not (rope straight): check distance 
end seesaw - end lid 

-> eventually set (rope straight) to TRUE 



move balloon 



output actual positions 



time loop's NEXT 



if max. time exceeded 



end 



Figure 7. Flow diagram for programming. 



5. Cases 

One further difficulty, which shall be called the "connecting-rod problem", has so far 
been overlooked. Of course, this will depend on the parameters, that is weights, 
lengths and initial angles. With the weight of the lid pulling down, one naturally 
expects the angular motion to slow down, that is yis and ifL become negative. The 
direction of in Fig. 6 because of the second contribution besides the one with 

(ps in it; '^TT^s^s'f's^Si will not necessarily be reversed. Furthermore, with directions 
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of forces at the end of the seesaw as in Fig. 6 there is no reason why we should 
not have > '^sars^os<ps _ xhen Fig.6 and (10) imply that < and 

^-^^^ > cosifL. (We shall not allow somersaults, that is we only regard ips < 
and (fL<7T/2. At (/^^ = 7r/2 the sandbag would drop from the seesaw anyway. And the 
lid of the box is assumed to engage when it reaches the vertical position. The projection 
factors sin a and sin(Q; + (fis ~ ^l) are always positive; the rope cannot come below 
the lid or above the right side of the seesaw.) A constellation with > ^gcosifL 

would mean that, instead of tearing, forces would be pushing into the rope from both 
ends. If then we solved (3) and (4) assuming the distance between the right end of the 
seesaw and the end of the lid constant, namely the length Ir, we would take the rope for 
a connecting rod and say that the lid having aquired enough angular momentum could 
then push against the seesaw from below via the "rope". In this case, in fact, the rope 
will loosen. The motion of the lid and the sandbag has to be separated, each merely 
further accelerated by the corresponding gravitational force (our massless seesaw giving 
way to the sandbag), until again the rope possibly straightens. This is why the general 
case of a horizontal velocity of the sandbag and the lid already moving had been treated 
at impact. The flow-diagram for integrating the equations of motion step by step in 
time in a numerical simulation is given in Fig. 7. 

6. Examples 

Four examples of calculations shall now be given. For each, parameters and starting 
values as defined in section 2 are given, dt is the time step for the calculation. The 
balloon is in contact with the hd of the box at the beginning. 

first example: Parameters: ms = 0.6 kg, th = 0.6 m, Iji = 1.0 m, tl = 0.2 m, 
mi = 0.5 kg, Xbox - Xseesaw = 0.73 m, hseesaw " hox = 0.33 m, dt = 0.001 s. 

Starting values: xs = -0.3 m, hs = 1.0 m, (ps = -25°, ^pi = 5°. 

Besides the horizontal and the vertical position of the sandbag in time, we have 
also plotted its trajectory in space. Seesaw and lid angles over time are united in one 
plot. The balloon's height as a function of time is drawn together with the height limit 
imposed by the hd. t—0, of course, is the time the sandbag starts falling. At t—0A2s it 
hits the seesaw. However, the rope is not straight yet, thus the seesaw just gives way 
to the falling sandbag. It already turns, such that at t=0.51s when the rope straigtens 
and the lid begins to be pulled up, the seesaw "starting" angle has changed from the 
initially given value. However, here, immediately after the lid is pulled up for the first 
time, the connection-rod check returns true, the rope loosens never to become straight 
again in this example. The most complicated differential equation (4) is only needed for 
a single time-step in this calculation! The one impact has given the hd enough angular 
momentum. It reaches the vertical position {(pL—7r/2) at i=0.82s. 
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t (s) t (s) 

Fig. 8. First example. 

Note that by this time the sandbag is on the seesaw's left arm 2 meters {^Jx^ + 
from the center. Our program does not contain a check for the sandbag to reach the end 
of the seesaw (which could be added without problems, nevertheless). In this example 
at later times the lid moves independently, and the sandbag dropping from the seesaw 
would not change that. The lid opens faster than the balloon can rise. Only shortly 
before the box has fully opened does the balloon touch the lid once again. With the 
lid in vertical position, the balloon can then rise freely. The time it takes to turn on 
the radio depends on the length of the balloon string and has not been calculated any 
more. Remark that the interesting part of the motion from the rope's first straightening 
until the lid's reaching vertical position only lasts from )f:=0.51s to t=0.82s, that is 0.31 
seconds. 

second example: As a second example a case has been sought where the rope stays 
straight at least for some time while pulling up the lid. Parameters and starting values 
are the same as in the first example except Ir = 0.62 m, hs = 0.5 m and ips = -35°. 
The sandbag meets the seesaw at t=0.24s, the lid starts moving at )f:=0.33s and reaches 
the vertical position at t=0.54s, that is already 0.21s later. The balloon cannot catch 
up with the lid during that time. At the end the sandbag is 0.6m from the center on the 
left arm of the seesaw, which is equal to the length of the right arm in this case. From 
t=0.49s on the sandbag cannot follow the seesaw any more and looses touch for more 
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than one time interval! (the rope loosens). We only show the two angles cps and cpL 
as functions of time, however, also compare them to the results that would have been 
obtained had the impact been calculated following conservation of the sum of angular 
momenta, energy or linear momenta (see section 3). These would be that the lid reaches 
the vertical position at t=0.50s, t=0.47s or t=0.34s, respectively. The last one definitely 
has to be disregarded as unphysical, though (see section 3). 




Fig. 9. Second example. 



third example: Parameters and starting values are the same as in the first example, 
except dt = 0.00025 s and at the beginning hs = 0.0 m, ips = 0°, that is the seesaw is 
horizontal and the sandbag is just placed on it, not let drop from a height. 




This example has been chosen, because the lid opens so slowly that the balloon 
catches up with it several times until practically staying trapped touching the lid. 
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Unfortunately, the chain reaction will not come to freeing the ballon and switching 
on the radio in this case. The seesaw will reach its vertical position at t=20s. However, 
there we assumed that the sandbag can slide an endless way on it, rs will be 1949m! The 
lid will be at v9l=74°. At (ps='n'/2 the sandbag finally dropping from the seesaw and the 
rope loosening, even if one persued this situation which has grown out of proportions, 
one constitutes that the hd does not have sufficient angular velocity to continue up to 
vertical position, but will quickly drop back down. 

forth example: Parameters: ms = 0.7 kg, th = 0.45 m, 1^ = 0.9 m, = 0.2 m, 

niL = 0.4 kg, Xbox - Xseesaw = 0.64 m, hseesaw " hox = 0.29 m, dt = 0.001 s. 

Starting values: xs = -0.3 m, hs = 0.5 m, ips = -20°, ipL = 5°. 

Whereas in the third example the rope 
stayed straight after the first impact 
at t=0.33s, here at t=0.40s the lid 
gets such a hard hit that immedi- 
ately afterwards the rope loosens again 
(connection-rod condition). The lid 
jumps up, falls back and only when 
its fall has again straightened the rope 
(t=0.88s) is pulled up by the counter- 
weight of the sandbag. However, as in 
the last example, the seesaw will reach 
vertical position at about t=10s (r^ will 
be 344m!), whereas the lid never will. 

The balloon catches up with the lid at t=0.64s, to be pressed down by it again and 
then follow it up again staying "glued" underneath it. 

7. Conclusions 

In conclusion Kermit's What-Happens-Next-Machine has been treated as an analytical 
mechanics and programming exercise. The two lever arms of the seesaw and the box's 
lid tied together by the rope can serve as a quite challenging example for finding the 
equation of motion using the Lagrange formalism. The most difficult point, however, is 
the decision which quantities are conserved and which are not at impact when the first 
force on the rope appears. A simulation has been set up and examples with different 
behaviour have been shown. It has been found that in cases where the counterweight 
of the falling sandbag is able to fully open the box to free the balloon, the interesting 
phase of the procedure only lasts some tenths of a second (this cannot be known from the 
original scene from Sesame Street, because there the machine does not work properly). 
Some simplifications have been made to be able to solve the problem theoretically (for 




0.0 0.5 1.0 1.5 2.0 



t(s) 

Fig. 11. Forth example. 
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example: no friction, neglection of the seesaw's mass). At the stage of development, 
including a mass for the seesaw would mean merely a small extension for who wants 
to refine the simulation. Attention has to be paid to the case that the seesaw then 
contributes inertia even if the sandbag is not in touch with it. As has been appreciated 
by some of our students, in the context of teaching theoretical mechanics, the What- 
Happens-Next-Machine can be given as an amusing tricky problem. 
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Attachments: 



exercise sheet from physics course (in German); skit of the scene is in English, 
program in FORTRAN, geometric angle and max. balloon height calculation 



References and Annotations 



[1] for example at NEMO in Amsterdam; 



http: / / video.aol.com/ video-detail/cliain-reaction-nemo-amsterdam/517940225 

[2] Domino Day for example. 

[3] just watch the scene yourself; http://www.youtube.com/watch?v=YnlrqVdFUOI| 
[4] Giiemez J, Fiolhais C and Fiolhais M 2009 Physics Education 44 53 

[5] In pratice the volume of helium has to be large enough to lift up the balloon with its skin which 
in turn requires a large enough box. The one in Fig.s 3 and 4 of about 25cmx25cm proved too 
small for that reason. 

[6] see any textbook on Theoretical or Analytical Mechanics, for example: 
H.Goldstein Klassische Mechanik, Aula-Verlag, Wiesbaden 1991. 

[7] see for example: Alonso/Finn Fundamental University Physics / chapter 6.4, Addison- Wesley 1980. 



Attachments 

Mutual dependence of the seesaw's and lid's angles and maximum height of the balloon 

Instead of the positions {xseesaw, hseesaw) and {xhox, h^ox) for the center of the seesaw 
and the hinge of the lid for our programs parameters / and 7 are chosen, that is the 
distance between these two points and the angle the line between them makes with the 
horizontal, rf, a, r and a are defined as can be seen from the drawing. Now the cosine- 
and sine-theorem can be used to get ifL when given ips- 
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d = \Jrjj + P - 2rHf cos{(ps + 7) , t = arcsm(rij sin((/?5 + 7) /d), 
cr = arccos((4r| + 0?^ — Z|j)/(4ri,o?)), </?i = 7r — 7 — r — cr, a = arccos((4r| + — 

rf')/(4rL/«)). 

The above calculation holds if the rope is straight. Having expressed d in terms of ips, 
the general constraint on (fi reads d"^ + (2r2,)^ — 2d ■ {2rL) ■ cos(7r — 7 — r — < /|.. 

As long as the lid has not come to vertical, the balloon must stay below it. We here 
calculate the maximum possible height of its center. There are shghtly different ways 
to do that, leading to the same result, of course. One could use the fact that the two 
drawn balloon radii together with the enclosed pieces of the box's wall and lid form a 
symmetric kite with two right angles. Here use the triangle with angle (pL formed to the 
left of the box by the extensions of the lid and the horizontal balloon radius. Putting 
together tarnpL = [hbox — hB^max) / and smipi = tbHtb + Ax) to eliminate Aa; one 
deduces + {hbox — hB,max)/ tampi = tb/ simpL and so the constraint on the balloon's 
motion is: hs < h^ox — ■ (1 — sinipi)/ cos ipL. 



What- Happens- Next- Machine 
Program 

program kermit 
implicit none 

integer cr , sos , steps , count , straight 

double precision hs, xs,phiS,phiL, gamma, f , pi, rL,rB,rH,ins,inL,iL,rs 
double precision rhoair ,rhohe , tmax ,dt ,g, vvs ,phisp,rsp , rspp, alpha 
double precision xf s ,hf s ,xf 1 ,hf 1 ,xsold, t ,philp, vhs 
double precision hb, vb,hbold,hbl ,hbf ,dphi ,phisold,philold 
double precision phispp,d,arccos,hsold, vvsold,vhsold,rsold,rspold 
double precision philplus.philminus, numerator, denominator 
double precision philstart, airline, lR,derivl,deriv2,phispold 
double precision pv,ph,coefa,coefb,coef c,coefd,philpold 
common pi, rH,rL,f, gamma, IR 

open ( 11 , f ile= ' whatnext . inp ' , status= ' old ' ) 
open (12 , f ile= ' whatnext . dat ' , status= ' unknown ' ) 

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



** This is what the input file ' whatnext . inp ' should look like ** 
************************************************************************ 

** * 

* 1.0 ; hs = initial height of sandbag in meters * 

* -0.3 ; xs = initial horizontal position of sandbag in meters * 

* -25.0 ; phiS = initial angle of the seesaw in degrees * 

* 5.0 ; phiL = initial angle of the lid in degrees * 

* 0.8 ; f = distance center of seesaw - lid's hinge in meters * 

* 1.0 ; IR = length of the rope in meters * 

* 27.0 ; gamma=angle of line seesaw center-lid hinge to horizontal (deg)* 

* 0.2 ; rL = half length of the lid in meters * 

* 0.6 ; rH = length of right side of seesaw in meters * 

* 0.6 ; ms = mass of the sandbag in kilograms * 

* 0.5 ; mL = mass of the lid in kilograms * 

* 1.29 ; rhoair = density of air in kilograms per cubic meter * 

* 0.18 ; rhohe = density of helium in kilograms per cubic meter * 

* 1.0 ; tmax = maximum time for calculation in second * 

* 0.001 ; dt = time step in seconds * 
** ** 
************************************************************************ 



readCll,*) 


hsold 


readCll,*) 


xsold 


readCll,*) 


phiS 


readCll,*) 


phiL 


readCll,*) 


f 


readCll,*) 


IR 


readCll,*) 


gamma 


readCll,*) 


rL 


readCll,*) 


rH 


readCll,*) 


ms 


readCll,*) 


mL 


readCll,*) 


rhoair 


readCll,*) 


rhohe 


readCll,*) 


tmax 


readCll,*) 


dt 


pi=4.0*atan(1.0) 


phiS=phiS/180.0*pi 


phiL=phiL/ 180 . 0*pi 


gamma=gamma/ 180 . 0*pi 


phisold=phiS 


philold=phiL 


philstart =phiL 


rB=rL 




iL=4 . 0/3 . 0*mL*rL*rL 


g=9.81 




vvs=0.0 




vvsold=0 . 




vhs=0 . 




vhsold=0 . 
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phisp=0 . 
philp=0.0 
phispold=0 . 
philpold=0.0 
vb=0 . 
phispp=0 . 
dphl=pi/200 . 

xf s=-f *cos (gamma) +rH*cos (phiS) 
hf s=f *s in (gamma) +rH* s in (phiS) 
xf 1=2 . 0*rL*cos (phiL) 
hf 1=2 . 0*rL*sin(phiL) 

airline=sqrt ( (xf s-xf 1) **2+ (hf s-hf 1) **2) 
if (airline .gt . IR) then 

print*, 'The rope is too short for these starting parameters. 

print* , airline 

stop 
endif 

if (phis. It. 0.0) then 

if (hsold. It . -xsold*tan(phis) ) then 

print*, 'sandbag must not start underneath seesaw' 

stop 

endif 
else 

if (hsold. It . xsold*taii(phis) ) then 
print* ,' sandbag must not start underneath seesaw' 
stop 
endif 
endif 
sos=0 

straight=0 
steps=int (tmax/dt) 
t=0.0 

if (abs(phil).le.l.OD-lO) then 

hbold=-rB 

else 

hbold=- (1 . 0-sin (phiL) ) /cos (phiL) *rB 

endif 

do count=l , steps 

t=t+dt 

cr=0 

c stopping condition 

if ((phiL. gt. pi/2. 0-dphi) .or. (phiS.gt. pi/2. 0-dphi)) goto 90 

c check if sandbag is on seesaw 
if (sos.eq.l) goto 20 

c sandbag and lid move separately 
vvs=vvsold+g*dt 
vhs=vhsold 

hs=hsold-vvsold*dt-0 . 5*g*dt*dt 
xs=xsold+vhsold*dt 
phiL=phiL+philp*dt 
philp=philp-mL*g*rL/iL*cos (phiL)*dt 
if (phiL. It .philstart) then 

phiL=philstart 

philp=0.0 
endif 

if (straight . eq. 0) then 
phiS=phiS+phisp*dt 
xf s=-f *cos (gamma) +rH*cos(phiS) 
hfs=f*sin (gamma) +rH*sin (phis) 
xf 1=2 . 0*rL*cos (phiL) 
hf 1=2 . 0*rL*sin(phiL) 

airline=sqrt ( (xf s-xf 1) **2+(hf s-hf 1) **2) 
if (airline. ge.lR) straight=l 
endif 

if (straight .eq. 1) then 
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call machphis(phiL,phiS) 
phisp= (phiS-phisold) /dt 
endif 

c check if sandbag has hit seesaw 

if (atan(hs/xs) .It.phiS) goto 30 

c initialize velocities at impact 
sos=l 

rs=-xs/cos (phiS) 

rsp=-vhs*cos (phiS) +vvs*sin (phiS) 
if (straight .eq. 1) then 
c keep angular momentum referred to ONE center 

d=sqrt (rH*rH+f *f -2 . 0*rH*f *cos (phiS+gamma) ) 

alpha=arccos ( (4 . 0*rL*rL+lR*lR-d*d) / (4 . 0*rL*lR) ) 

call machphiLCphiS ,phiL) 

call machphiL(phiS+dphi .phiLplus) 

call machphiL(phiS-dphi ,phiLniinus) 

derivl= (phiLplus-phiLminus) /2 . 0/dphi 

numerator=ms*rs* (vvs*cos (phiS)+vhs*sin(phiS) )+(iL*philp*rH* 
1 sin (alpha+phiS-phiL) ) / (2 . 0*rL*sin (alpha) ) 

denominator=ms*rS*rS+(iL*rH*sin(alpha+phiS-phiL)*derivl)/ 
1 (2.0*rL*sin(alpha)) 

phi sp=numer at or/ denominator 

philp=derivl*phisp 
cc if we kept sum of angular momenta 

cc sandbag around center of seesaw and lid around its hinge 



c call machphiL(phiS,phiL) 

c call machphiL(phiS+dphi,phiLplus) 

c call machphlL(phiS-dphi,phlLminus) 

c derivl= (phiLplus-phiLminus) /2 . 0/dphi 

c phisp=(ms*rS*(vvs*cos(phiS)+vhs*sin(phiS))+iL*philp)/ 

c 1 (ms*rS*rS+iL*derivl) 

c philp=derivl*phisp 

cc if we kept kinetic energy 

c call machphiL(phiS,phiL) 

c call machphiL(phiS+dphi .philplus) 

c call machphiL(phiS-dphi ,philminus) 

c der ivl= (philplus-philminus ) /2 . 0/ dphi 

c nuinerator=ms* (vvs*cos (phiS) +vhs*sin(phiS) ) **2+iL*philp*philp 

c denominator=ms*rS*rS+iL*derivi*derivl 
c phisp=sqrt (numerator/denominator) 

c philp=derivl*phisp 

cc if we kept linear momentum - results in nonsense 

c ph=ms*vhs-mL*rL*philp*sin(phiL) 

c pv=-ms*vvs+mL*rL*philp*cos (phiL) 

c call machphiL(phiS,phiL) 

c call machphiL(phiS+dphi, philplus) 

c call machphiL(phiS-dphi,philminus) 

c derivl= (philplus-philminus) /2 . 0/dphi 

c coef a=-ms*cos(phiS) 

c coef b=ms*rS*sin (phis) -mL*rL*sin(phiL)*derivl 

c coef c=-ms*sin (phis) 

c coef d=-ms*rS*cos (phiS)+mL*rL*cos(phiL)*derivl 

c rsp=(ph*coef d-pv*coefb)/(coef a*coef d-coefb*coef c) 

c phisp= (ph*coef c-pv*coef a) / (coef b*coef c-coef a*coef d) 

c philp=derivl*phisp 
else 

phisp=vvs*cos (phiS) /rs+vhs*sin(phiS) /rs 
endif 



c anyway, that's it for this time step 
goto 30 

20 continue 
c sandbag is on seesaw 

if (straight . eq. 1) then 

call machphiL (phiS , phiL) 

call machphiL (phiS+dphi .philplus) 

call machphiL (phlS-dphi ,philminus) 

derivl= (philplus-philminus) /2 . 0/dphi 
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deriv2=(philplus-2 . 0*phiL+philinlnus) /dphi/dphi 
numerator=ms*g*rS*cos (phiS) -2 . 0*ms*rS*rsp*phisp- 
1 mL*g*rL*cos (phiL) *derivl-iL*phisp*phisp*deriv2*derivl 
denominator=ms*rS*rS+iL*derivl*derivl 
phi spp=numerator /denominator 
else 

phispp= (g*cos (phiS) -2 . 0*rsp*phisp) /rS 
endif 

rspp=rs*phisp*phisp+g*sin(phiS) 

rsp=rsp+rspp*dt 

rS=rS+rsp*dt 

phisp=phisp+phispp*dt 

philp=derivl*phisp 

phiS=phiS+phisp*dt 

xs=-rS*cos (phiS) 

hs=-rS*sin(phiS) 

vvs=rsp*sin(phiS) +rS*cos (phiS) *phisp 
vhs=-rsp*cos (phiS) +rS*sin(phiS) *phisp 
if (straight . eq. 1) then 
call machphiL (phiS , phiL) 
else 

phiLp=philp-mL*g*rL/iL*cos (phiL)*dt 

phiL=phiL+philp*dt 

endif 

if (phiL. It .philstart) then 

phiL=philstart 

philp=0.0 
endif 

check if sandbag falls faster than allowed by free fall 
if (vvs-vvsold.gt .g*dt) then 
sos=0 

straight=0 
t=t-dt 
hs=hsold 
xs=xsold 
vhs=vhsold 
vvs=vvsold 
rs=rsold 
rsp=rspold 
phiS=phisold 
phiL=philold 
phisp=phispold 
philp=philpold 
goto 100 
endif 

if (straight. eq.l) then 
connection rod check, check whether rope loosens 
if (ms*rS*rS*phispp+2 . 0*ms*rs*rsp*phisp . gt . 
1 ms*g*rS*cos(phiS)) then 
cr=l 

straight=0 
endif 
else 

distance check, check whether rope straightens 
print* ,' check distance' 
xfs=-f*cos (gamma) +rH*cos (phiS) 
hf s=f *sin (gamma) +rH*sin (phiS) 
xf 1=2 . 0*rL*cos (phiL) 
hf 1=2 . 0*rL*sin(phiL) 

airline=sqrt((xf s-xf l)**2+(hf s-hfl)**2) 
if (airline .ge . IR) then 
straight=l 
endif 
endif 

30 continue 
move the balloon 

hbl=- ( 1 . O-sin(phiL) ) /cos (phiL) *rB 
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hbf =hbold+vb*dt+0 . 5*g* (rhoair-rhohe) /rhoair*dt*dt 

if (hbl.lt.hbf ) then 

hb=hbl 

vb=0 . 

else 

hb=hbf 

vb=(hb-hbold)/dt 
endif 

c update values 
hbold=hb 
xsold=xs 
hsold=hs 
vvsold=vvs 
vhsold=vhs 
rsold=rs 
rspold=rsp 
phisold=phiS 
philold=phiL 
phispold=phisp 
philpold=philp 

c output 

write (12 ,123) t ,xs ,hs ,phiS,phiL,hbl,hb,dble(cr) ,dble( straight) 
123 format (9F8.3) 

100 enddo 

print*, 'time out' 
goto 95 

90 print*, 'reached vertical position, t',t 

95 close(ll) 
close(12) 
end 

double precision function arccos(x) 
implicit none 

double precision x,pi,tangens 

pi=4.0*atan(1.0) 

if (abs(x) .le.l.OD-10) then 

arccos=pi/2.0 
else 

tangens=sqrt ( 1 . 0-x*x) /x 

arccos=atan(tangens) 

if (x. It. 0.0) arccos=arccos+pi 
endif 
return 
end 

subroutine machphiS(phiL,phiS) 
implicit none 

double precision phiS, phiL, pi, lR,rH,rL,f, gamma, dhere,tausigma 
double precision arccos ,downangle ,upangle 
common pi ,rH,rL , f ,gamma,lR 
tausigma=pi-phiL-gamma 

dhere=sqrt (f *f +4 . 0*rL*rL-4 . 0*f *rL*cos (tausigma) ) 

downangle=arccos ( (dhere*dhere+f *f -4 . 0*rL*rL) /2 . 0/dhere/f ) 

upangle=arccos ( (dhere*dhere+rH*rH-lR*lR) /2 . 0/dhere/rH) 

phi S=do wnangl e+upangl e -gamma 

return 

end 

subroutine machphiL(phiS,phiL) 
implicit none 

double precision pi ,rH,rL, gamma, lR,f ,phiL,phiS, arccos, tau, sigma,d 

common pi, rH,rL,f, gamma, IR 

d=sqrt (rH*rH+f *f -2 . 0*rH*f *cos (phiS+gamma) ) 

tau=arccos ( (f *f +d*d-rH*rH) /2 . 0/f /d) 



What- Happens- Next- Machine 



20 



sigma=arccos ( (4 . 0*rL*rL+d*d-lR*lR) /4 . 0/rL/d) 

phiL=pi-gamma-tau-sigma 

return 

end 



7.1. Skit of the scene from Sesame Street 

(Kermit is standing neoct to a rope on a pulley that's attached to a sandbag. The sandbag 
is suspended over a seesaw. The other end of the seesaw, pointing down, has a string 
on it attached to the closed lid of a box containing an inflated helium balloon. Its string 
is trailing out of the box and is tied to the switch of a radio.) 

Kermit (who has a pair of scissors): First of all, I'll cut the rope. Then the sandbag 
will fall on the seesaw. Then that end of the seesaw will go down, and the other end of 
the seesaw will go up. ...and the balloon will float up in the air! And you will notice 
that the balloon is tied to the switch of the radio. 

Kermit (going over to the radio): You see, in the old days I had to walk all the way 
over here to turn on my radio. It would tire my flippers. It would waste time. But not 
anymore ... thanks to the magic of What-Happens-Next! ... First I cut the rope! (He 
does ... and the sandbag stays suspended in mid-air.) 

Kermit: And what happens next is ... the sandbag seems to be stuck! The rope 
must be tangled ... (He gives it a few good whacks, but it still doesn't fall.) Oh, well, 
we'll just skip the sandbag part and go on to the seesaw. I'll just give it a push. (He 
does, but it won't budge.) The seesaw seems to be stuck ... (He gives a few determined 
grunts and shoves, but to no avail.) 

Kermit: Oh, well, we'll just ignore the stuck seesaw and move to the box. And 
what do you know, I'll bet the lid won't open. (He gets it open on the first try, though.) 
Oh, there. (The balloon refuses to rise.) I'll give it a few good kicks! (He does, but the 
balloon stays put. By now he 's really upset.) Come on, balloon! Are you going to float 
up in the air or not? (Kicks it again.) All right, balloon, I'm giving you one last chance! 
You float up in that air right now, or 1 won't be responsible for what happens next! (He 
gives the box a vicious kick, and finally the balloon floats up in the air, turning on the 
radio, which starts playing music.) 

Kermit: Thanks to the magic of What-Happens-Next. Once I get all the bugs out. 
(Just then, he sees the balloon carrying the radio away.) 

Kermit: What happens next is I gotta get another radio! (He jumps as the sandbag 
suddenly falls on the seesaw, which breaks in two.) 

[adapted from SilveryShoe@aol.com] 



